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Abstract 

In this work, using the scattering matrix method, we have investigated the transmis- 
sion coefficients and the thermal conductivity in a double-bend waveguide structure. 
The transmission coefficients show strong resonances due to the scattering in the 
midsection of a double-bend structure; the positions and the widths of the reso- 
nance peaks are determined by the dimensions of the midsection of the structure. 
And the scattering in the double-bend structure makes the thermal conductivity 
decreases with the increasing of the temperature first, then increases after reaches 
a minimum. Furthermore, the investigations of the multiple double-bend structures 
indicate that the first additional double-bend structure suppresses the transmission 
coefficient and the frequency gap formed; and the additional double-bend structures 
determine the numbers of the resonance peaks at the frequency just above the gap 
region. These results could be useful for the design of phonon devices. 
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1 Introduction 



Since the discovery of the quantized electronic conductance phenomena [1,2], 
the double-bend electron waveguide has been investigated as a very impor- 
tant case. The problem was firstly studied by Weisshaar et al. [3] using mode- 
matching technique. The experiment work on the low temperature conduc- 
tance of the double-bend waveguide was carried out by Wu et al. [4] By 
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using the recursive Greens function technique, Kawamura et aL [5,6] stud- 
ied this double-bend waveguide again. All the works showed strong resonant 
transmission due to internal reflections in the special structure and Refs. [5,6] 
showed the "existence of an energy gap between the first and second subband 
threshold energies where the conductance is suppressed for multiple double- 
bend structures" . Based on the property of the strong resonant transmission 
through the double-bend waveguide, Shi et al. [7] proposed a simple spin filter 
with which an extremely large spin current is expected. 

Same as the electronic conductance, the thermal conductivity is also important 
for semiconductor nanostructures. For an ideal elastic beam at an enough low 
temperature, the thermal conductivity is dominated by the ballistic phonon 
and is quantized in a universal unit, n'^k%T/3h, analogous to the well-known 
2e^//i electronic conductance quantum. [8,9,10] These predictions have been 
verified experimentally [11]. Since then, many works have been done to study 
the geometrical effects on the transmission of the acoustic phonon and the 
thermal conductivity through the quantum waveguide with various geometries 
[12,13,14,15,16,17,18,19,20]. But as an important case, the quantum waveguide 
with a double-bend structure has not been studied on the phonon transmission 
coefficients and thermal conductivity. Therefore, in this work, we calculate the 
phonon transmission coefficients and thermal conductivity in a double-bend 
quantum waveguide using the scattering-matrix method [21,22,23,24,25] and 
considering the stress-free boundary condition [15,16,17,18,19,20]. 

The organization of this paper is as follows. In section 2, we present the model 
and the numerical method briefly. The numerical results for a one double-bend 
structure are discussed in section 3. In section 4, the results of a two double- 
bend structures are presented. Finally, a summary is made in section 5. 



2 Model and formalism 

The geometry of the double-bend quantum waveguide is sketched in flgure 1. 
Regions-I and -III are the leads of the device; region-II is the midsection. 
The parameters of this device are W, L and Di, which are the lateral width 
of the leads, the lateral width of the midsection and the longitudinal length 
of the midsection, respectively. We assume that the temperature in the two 
leads (regions-I and -III) are Ti and T2; and the temperature difference 6T 
(ST — Tl — T2 > 0) between two leads is very small. The mean temperature 
T {— {T1+T2) /2) may then be adopted as the temperature of two leads. In this 
work, we choose the same thickness for the three regions and let the thickness 
smaller than the other two dimensions and also than the wavelength or the 
coherence length of the elastic waves; there is no mixing of the z modes and 
a two-dimensional calculation is then adequate [8,9,10,15]. For the imperfect 
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Fig. 1. (a) Double-bend waveguide with lead width W, lateral length L, and longi- 
tudinal length Di . (b) Two double-bend waveguides in a series with the longitudinal 
length D2. 

contact at the regions-I and -III, the thermal conductance K at temperature 
T is given by [8,15] 

K = - — =7: — / Tmi'^) rjda;, (1) 

where {u) is the transmission coefficient from mode m of region-I at fre- 
quency uj across all the interfaces into the modes of region-III; uj^ is the cutoff 

frequency of the mth mode; (3 = 1/ (k^T), is the Boltzmann constant, T 
is the temperature; and h is the Plank's constant. It can be seen from Eq. (1) 
that the central issue of the problem is to obtain the transmission coefficient, 

Wc consider the scalar model for the clastic wave; the model for thin geometry 
at low temperature is used so that the calculation is two-dimensional. Here, 
we treat the simplest case for a horizontally polarized wave (SH) (polarized 
along the z direction) propagating in the x-direction. So the wave equation of 
the displacement field u is 

^ - vl,V\ = 0, (2) 
where vsr is related to the mass density p and elastic stiffness constant C44 

by 

^^SH = \/c^47p- (3) 
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The stress-free boundary conditions at the edges require n ■ Vu = 0, where n 
is the unit vector normal to the edge. For the double-bend structure depicted 
in Fig. 1, the solution to Eq. (2) in region-,^ (regions I, II and III) can be 
expressed as 



N 



m=0 



0m (y) , 



(4) 



where is the reference coordinate along the x-direction in region-,^; k^^ is 
the wavenumber of the transmitted and reflected waves in region-,^, giving by 
the energy conservation condition 

u;'^km^' + m\\l^'/Wl, (5) 

with the transverse dimension of region-^; 0^ (y) represents the transverse 
wavefunction of acoustic mode-m in region-^, 



<I>1 



(m = 0) ; 

2 cos ^y (m ^ 0) , 
(m = 0) ; 



^cos^{y-L) (m^O), 



(6) 
(7) 



In principle, the sum over m in Eq. (4) includes all propagating modes and 
evanescent modes (imaginary k^). However, in the practical calculation, be- 
sides all the propagating modes we only take a same limited number of evanes- 
cent modes into account to meet the desired precision for each region in a 
double-bend structure. The boundary matching conditions require the con- 
tinuity of the displacement u and the stress C/^^du/dx at the interface of 
regions-I and -II and the interface of regions-II and -III. So we can obtain 
the equations for the coefficients in Eq. (4). Rewriting the resulted equations 
in the form of matrix, we can derive the transmission coefficient, r^, by the 
scattering matrix method [21,22,23,24,25]. 

In the calculations, we employ the following values of elastic stiffness constant 
and the mass density [26]: C44(GaAs) = 5.99(10^°Nm-2) and p(GaAs) = 
5317.6 (kgm~^), and choosing 1^ = 10 nm. 
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3 Numerical results for a one double-bend structure 




(b) 

Fig. 2. Transmission coefficient versus the lateral length L and the incident 

phonon frequency w/Agu for the double-bend structure with W = 10 nm, where 
^SH = f^m+i —^m = T^vsa/W (vsu IS the acoustic wave velocity in GaAs) represents 
the splitting of the cutoff frequency between the (m + l)th and mth mode in the 
region-I. (a) Di = 10 nm; (b) Di = 12 nm. 



For acoustic phonon mode, the stress-free boundary condition allows acous- 
tic waves propagate through the structure in the zero-th mode with cutoff 
frequency uj = 0. First, we discuss the lowest SH acoustic mode propagate 
through the double-bend structure shown in Fig. 1(a) from region-I to region- 
Ill. The transmission coefficients of the zero-th mode versus the lateral length 
L and the incident phonon frequency a;/ Ash for the double-bend structure 
with W = 10 nm are shown in Figs. 2(a) {Di — 10 nm) and 2(b) {Di — 12 nm). 
One can find from Fig. 2 that the transmission coefficient approaches unity as 
u ^ and is independent of the lateral length L. This is because, at u; ^ 0, 
the wavelength of the phonon is much larger than the dimension of the lateral 
length L of region-II and the displacement field u then becomes essentially 
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the same throughout. Wc can also find from Fig. 2 that, for L = 10 nm, the 
transmission coefficient is unity for any incident frequency. In this case the 
double-bend waveguide reverts to be a uniform waveguide; and in the uniform 
waveguide, the transmission coefficient of the zero-th mode of phonon is unity. 
Furthermore, from Fig. 2 one can find more resonance like features occur in 
the transmission coefficient with increasing incident frequency or the lateral 
length of region-II, L; and the influence of the longitudinal width Di also 
appear. 

To clearly see the resonance structure in the transmission coefficient, the trans- 
mission coefficient as a function of cj/ Ash for different lateral lengths of region- 
II, L, and as a function of L for different a;/ Ash are shown in Fig. 3. 

In region-II, the propagation wave and reflected wave with the same wavenum- 

ber, k — ^Jlo"^ /v'^ — m^TT^/L*!^, will couple together. The reflected wave from 
region-II will interfere with the propagation wave in region-I. When the phase 
shift between the incident wave in the region-I and the reflected wave from 
the region-II reaches 2n7r, the transmission will be enforced, and, tq = 1. For 
a flxed lateral length L, same as the resonant transmission of electrons [27], 
the resonance condition is approximately 2A;„L = 2n7r. In Figs. 3(a) and 3(b), 
the arrows mark the resonance positions. In Fig. 3(b), when a;/Agjj < 1, the 
resonance peak positions are similar as the positions in Fig. 3(a). This is easy 
to understand: in this case only the zero-th mode can exist in region-II, so 
the wavenumber k = lu/v is irrelevant with Di. But for a;/Agjj > 1, because 
of more than one transverse mode existing in the region-II the transmission 
coefficients in Fig. 3(b) are totally different with those in Fig. 3(a). Further- 
more, from Figs. 3(a) and 3(b), one can find that the closer the frequency to 
the cutoff frequency of the 1th mode in region-II, the narrower the resonance 
peak is. This is due to the fact that at the cutoff frequency of the 1th mode, 
the 1th mode is excited and the resonant condition is destroyed totally. The 
transmission coefficients as a function of the lateral length L show periodic 
resonance in Figs. 3(c) and 3(d). The period is determined by the condition 
2k AL = 2tt as, AL = ir/k — X/2. These results could be useful for the design 
of phonon devices. 

From the transmission coefficients the thermal conductivity can be obtained 
from Eq. (1). Let u; — rj Ask, T — CAsh^/^b and choose h—1'm Eq. (1), then 



where Gq = ir'^k'^T/Sh is the universal value of the thermal conductivity. Let's 
consider the integration in above equation: if m = 10 and C ^ 1) then rj^ = 10, 
eVC > e^o > 1, and 
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Fig. 3. Transmission coefficient as a function of: (a) the reduced frequency cj/Ash 
with W = 10 nm and Di = 10 nm; (b) the reduced frequency w/Agn and w/Agjj 
with W = 10 nm and Di = 12 nm, where u/Ag^ = -kvsk/Di. Curves are vertically 
offset by 1.5 units for clarity. The arrows indicate the resonance peaks. Transmis- 
sion coefficient as a function of the lateral length L with: (c) Di = 10 nm; (d) 
Di = 12 nm. The dotted, dashed and solid curves in (c) and (d) correspond to 
a;/AsH = 0.2, 0.5 and 0.8, respectively; and in (d), these curves correspond to 
oj/A^^ = 0.24, 0.6 and 0.96, respectively. 
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in which we have used the equation, 



f 

Jo 



e ^'^x^'dx = 



,n+l ■ 



n = 0, 1,2, s>0. 



Therefore, the sum of m in Eq. (1) can truncated at m = 9 for ^ < 1. 
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Fig. 4. Thermal conductivity reduced by G^, as a function of the reduced temper- 
ature, C = ^B^/^^SHi for the different lateral lengths of L: (a) the total thermal 
conductivity; (b), (c) and (d) are the thermal conductivity of the modes-0, 1 and 2, 

respectively. 

Fig. 4 shows the reduced thermal conductivity K/Gq as a function of the 
reduced temperature k-^T / h/S.^]i- One can find from Fig. 4(a) that the to- 
tal reduced thermal conductivity approaches to the unity (i.e. the thermal 
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Fig. 5. Transmission coefficient as a function of the reduced frequency a;/ Ash- (a), 
(b) and (c) are for tq, ti and T2, respectively. 

conductivity, K/T, approaches to the ideal universal value, n'^k'^/Sh) when 
T — > 0; the result is independent of the lateral length of the double-bend 
waveguide. This is because the wavelength of the incident phonons is much 
larger than the dimension of the region-II when T — * and, thus, the scat- 
tering of phonons in the region-II is negligible. Fig. 4(a) shows a thermal 
conductivity plateau at L = and at low temperatures due to no scattering; 
and the thermal conductivity increases lineally with the temperature because 
the higher transverse modes are excited at higher temperatures (as shown in 
Figs. 4(c) and 4(d)). From Fig. 4(a), for L 7^ 0, one can find the thermal 
conductivity decreases from the universal value to a minimum at low temper- 
atures, then increases with the the temperature. This can be explained from 
Figs. 4(b)-4(d): With increasing temperature, the thermal conductivity of the 
mode-0 decreases for L 7^ 0, so the total thermal conductivity decreases from 
the universal value. When the temperature further increases, the higher mode 
can be excited and the total thermal will increase after reach the minimum. 
Furthermore, Fig. 4(a) shows that the total thermal conductivity decreases 
with increasing the lateral length; and when L > 50 nm, the total thermal 
conductivity is insensitive to further increasing of L. 

From Fig. 5(a), one can find that the transmission coefficient of the zero-th 
mode. To, for L — 20 nm is higher than that for L — 50 and 90 nm; and tq 
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for L = 50 nm is similar as for L = 90 nm. So for the thermal conductivity 
of mode-0 shown in Fig. 4(b), it is higher for L = 20 nm than for L = 50 
and 90 nm. Fig. 5(b) shows the transmission coefficients of the 1th mode for 
different lateral lengths of L. Here, Ti for L = 90 nm is higher than that for 
L = 50 nm, so the thermal conductivity for L = 50 nm is lower than that for 
L = 90 nm in Fig. 4(c). The thermal conductivity for L = 50 nm is higher than 
that for L = 20 nm shown in Fig. 4(d), because the transmission coefficient 
of 2th mode for L = 50 nm is higher than that for L = 20 nm in Fig. 5(c). 



4 Numerical results for two double-bend structures 
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Fig. 6. (a), (b) Transmission coefficient as a function of the reduced frequency a;/ Ash 
for W = 10 nm, L = 20 nm and Di = 10 nm. (a) A two double-bend structure for 

different longitudinal lengths of D2, where curves are vertically offset by 1.5 units 
for clarity, (b) A two double-bend structure and a three double-bend structure for 
the same value of D2. 



Fig. 6 shows the transmission coefficient through a quantum waveguide with 
more than one double-bend structures. Fig. 6(a) represents the case of a two 
double-bend waveguide for different longitudinal lengths of One can find 
from the figure that an additional double-bend structure makes the frequency 
gap appear; and the gap width decreases with D2 increases. Furthermore, it is 
obviously that the width of resonance peak at the frequency just above the gap 
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region decreases with D2 increases. Fig. 6(b) compares the transmission coef- 
ficients through a two double-bend structures and a three double-bend struc- 
tures for the same value of Adding an additional double-bend structure to 
the two double-bend structures doesn't change the frequency gap width but 
makes an additional resonance peak appearing at the frequency above the gap 
region. It may draw conclusions that the first additional double-bend structure 
to the single double-bend structure suppresses the transmission coefficient and 
forms a frequency gap; and an additional resonance peak is formed at the fre- 
quency just above the gap region for each additional double-bend structures. 
We can easily expect that an additional double-bend structure increases the 
scattering of the phonons so that the thermal conductivity decreases. 



5 Summary 

Using the scattering matrix method, we have investigated the phonon trans- 
mission and thermal conductivity in a double-bend waveguide structure. It is 
observed strong resonant transmission determined by the scattering in region- 
II. The position and the width of the resonance peaks are determined by the 
longitudinal length Di of region-II. The thermal conductivity approaches to 
universe value for T ^ 0. This is because the wavelength of the incident 
phonons is much larger than the dimension of the region-II when T ^ 0; 
the scattering of phonons by the region-II is negligible. Due to the scatter- 
ing in region-II, the thermal conductivity of mode-0 decreases with increasing 
temperature; and at high temperatures, more than one mode are exited and 
the thermal conductivity of these exited modes increases. The total thermal 
conductivity decreases thus to a minimum, then increases with the tempera- 
ture. The transmission coefficient of a double-bend structure determines the 
thermal conductivity for different lateral lengths of L. 

In this paper, the transmission coefficients and the thermal conductivity in 
multiple double-bend waveguide structures are also studied. The results show 
that the first additional double-bend structure to the single double-bend struc- 
ture suppresses the transmission coefficient and forms a frequency gap; and 
an additional resonance peak appears at the frequency just above the gap 
region for each additional double-bend structure. The additional double-bend 
structure suppresses more the thermal conductivity. 
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